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Abstract. Structure of weakly bound/unbound nuclei close to particle drip lines is different 
from that around the valley of beta stability. A comprehensive description of these systems 
goes beyond standard Shell Model (SM) and demands an open quantum system description 
of the nuclear many-body system. We approach this problem using the Gamow Shell Model 
(GSM) which provides a fully microscopic description of bound and unbound nuclear states, 
nuclear decays, and reactions. We present in this paper the first application of the GSM for a 
description of the elastic and inelastic scattering of protons on 6 He. 



1. Introduction 

The nuclear SM was proposed almost sixty years ago[H [2J. Soon after, the interacting SM 
was developed and used extensively to understand a wealth of data on nuclear levels, moments, 
collective excitations, electromagnetic and j3 decays, and various particle decays. Interacting SM 
describes the nucleus as a closed quantum system: nucleons occupy bound, hence well localized, 
single-particle orbits of a harmonic oscillator potential and are isolated from the environment of 
unbound scattering states. Since the scattering continuum is not considered explicitly, the 
presence of decay thresholds and the interplay of Hermitian (internal) and anti-Hermitian 
(external, via the continuum) mixing of SM configurations is neglected. This competition 
yields a complicate interference pattern[3] and is a source of many collective features such as, 
e.g., the resonance trappingJH [5j [6] and super-radiance phenomenon^ [8], the multichannel 
coupling effects in reaction cross-sections [9 J and shell occupancies [10j. the clustering[llj, the 
modification of spectral fluctuations [12], the deviations from Porter-Thomas resonance widths 
distribution [6, 13J, and so on. 

It was clear in fifties that the SM is a phenomenological tool which employs the strongly 
renormalized bare nuclear interaction between nucleons |14j and neglects coupling to the 
continuum [15]. The overwhelming success of the SM, its elegance and simplicity resulted in 
neglecting most of these questionable assumptions. Unfortunately, the divide between discrete 
and continuum states in SM has led to an artificial separation of nuclear structure from nuclear 
reactions, and hindered a deeper understanding of nuclear properties. It is often believed that 
the understanding why the simple interacting SM works so well will be advanced when the goal 
of ab initio many-body approaches including the continuum coupling will be achieved. 



Many structural properties of the nucleus are determined by means of nuclear collisions and 
this calls for a unified theoretical framework. Feshbach at the beginning of sixties formulated a 
unified theory of nuclear reactions using the effective Hamiltonian and the projection operator 
method to select the open channel components of the wave function[16j. This development led 
to various formulations of the real-energy continuum shell model (CSM) [T71 118 } [T9l I20|, |2T| I22j . 
Nowadays, the real-energy CSM provides a unified description of the structure and reactions 
with up to two nucleons in the scattering continuum[2Ul |2"2"] . 

At the same time, Fano noticed that the exact coincidence of different configurations above 
the lowest particle emission threshold makes the perturbation theory inadequate and calls for a 
generalization of the standard SM[23J. The achievement of Fano's goal took almost forty years 
and required the development of new mathematical approach of the rigged Hilbert space [24J, 
new methods to deal with diverging integrals (matrix elements of one- and two-body operators) 
involving resonance and scattering states [281 [291 [301 EH [32], the formulation of the generalized 
completeness relation including single-particle (s.p.) bound states, resonances and scattering 
states[25j. These different and independent developments in mathematics and physics enabled 
finally a satisfactory formulation of the new many-body theory, the GSM[33l [Ml (33 [36], which 
offers a fully symmetric treatment of bound, resonance and scattering states in the multiparticle 
framework. In this formulation, the maximum number of particles in the scattering continuum 
is not a priori prescribed as in the real-energy CSM, but follows from the Schrodinger variational 
principle for the many-body Hamiltonian. 

Until now, GSM has been primarily used in the context of nuclear structure. (For a recent 
review, see Ref. [26].) In this paper, we shall extend GSM to reaction problems by the coupled- 
channel (CC) formulation of the scattering problem. The very first application of the GSM-CC 
formalism will be presented in this paper for the proton scattering on 6 He target. The proposed 
GSM-CC formalism can be easily straightforwardly generalized for the description of nuclear 
reactions in the ab initio framework of the No-Core Gamow Shell Model [27], 

In Sect. 2 we present essential features of the GSM. The reaction wave function and the 
derivation of GSM-CC equations in the coordinate space representation are discussed in Sect. 
3. The method for solving the CC equations is presented in Sect. 4. The first application of the 
GSM-CC formalism for a description of 6 He+p reaction is discussed in Sect. 5. Finally, main 
conclusions of this work are summarized in Sect. 6. 

2. The Gamow Shell Model 

Resonance phenomena are generic intrinsic properties of many mesoscopic systems. What is 
specific to atomic nuclei are strong nucleon-nucleon correlations which impose a simultaneous 
description of the configuration mixing and the coupling to decay channels. The fundamental 
difficulty in the CSM formulation and the reason why the Fano's program [23J had to wait forty 
years to find a comprehensive solution is the fact that resonances do not belong to the Hilbert 
space. Therefore, whatever strategy is chosen to formulate a configuration interaction approach 
in open quantum systems, a key points are always: (i) the treatment of s.p. resonances in a 
many-body framework, and (ii) the definition of complete s.p. and many-body bases. 

For the GSM, a s.p. basis is given by the Berggren ensemble [25] which consists of Gamow 
(resonant) states and the non-resonant continuum (see Fig. [1]). (For a detailed description of 
the GSM see Ref. [26].) The GSM Hamiltonian is Hermitian. However, since the s.p. vectors 
have either outgoing or scattering asymptotics, the Hamiltonian matrix in GSM is complex 
symmetric and its eigenvalues are complex above the first particle emission threshold. Hence, 
both real-energy and complex-energy CSM formulations lead to a non-Hermitian eigenvalue 
problem above the threshold and contain all salient features of an interplay between opposite 
effects of Hermitian and anti-Hermitian couplings. 
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Figure 1. The schematic illustration of a complete s.p. basis in fc-plane for GSM calculations. 
The complete many-body basis is spanned by Slater determinants involving nucleons on bound, 
resonant and continuum s.p. states along the contour L+ in the complex fc-plane. GSM 
Hamiltonian matrix is complex symmetric and yields complex energy eigenvalues above the 
first particle emission threshold. 



2.1. The Hamiltonian 

The translationally invariant GSM Hamiltonian in intrinsic nucleon-core coordinates of the 
cluster-orbital shell model|37|, can be written as: 
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where M c is the mass of the core, Hi is the reduced mass of either the proton or neutron 
(1/fjLi = 1/raj + 1/M C ), U is the s.p. potential describing the field of the core, V is the two-body 
residual interaction between valence nucleons. The last term in (pQ) represents the recoil term. 

The particle-core interaction is a sum of nuclear and Coulomb terms: U = U + U . 
The nuclear potential U is approximated either by a Woods-Saxon (WS) field with a spin- 
orbit term[35j or by the Gamow-Hartree-Fock (GHF) potential [36J. The Coulomb field U is 
generated by a Gaussian density of Z c core protons [38]. Similarly, the residual interaction can 
split into nuclear and Coulomb parts: V = V N +V ' , where V N is the modified surface Gaussian 
interaction (MSGI)[38| and V c is the two-body Coulomb interaction. 
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asymptotic behavior of the Coulomb interaction. The second term in this equation and the 
two-body recoil term can be expanded in the harmonic oscillator basis [39, 38^ which provides an 



accurate treatment of the long-range physics of the Coulomb potential. In this work, we took 9 
harmonic oscillator shells with the oscillator length 6 = 2 fm. 

3. N-body GSM reaction wave functions 

The CC framework is a convenient to formulate the GSM description of reactions involving 
one proton (neutron) scattering processes. The CC equations are obtained from the following 
A-body matrix elements: 

{A{(* f \ J ' ® (r f l f j f T f \}j$ A \H\Am) Ji ® \n i { ji n)}{j A ) , (2) 

where (I^/)m ) are m itial (final) GSM eigenvectors of (A — l)-body system, r*, rf are 

radial coordinates, r,, Tf are isospin quantum numbers (proton or neutron), and £%, ji, £f, jf are 
angular quantum numbers. All ^4-body wave functions are fully antisymmetrized, as emphasized 
by the A symbol. 

In order to express the antisymmetry in a convenient way, the \r£jr) channel is expanded in 
a s.p. basis of GSM wave functions u n (r) generated by the s.p. potential Ubasis- This implies 
for the associated creation operator: 

a Ujr = J2 u n(r) a{ ljT . (3) 

n 

Hence, the expression of considered ^4-body wave function becomes: 

A\m J ®\rljr)\^ A ) = £u„(r){a^ r |*> J }^ . (4) 

n 

States of the target [A — l)-nucleus are eigenstates of the Hamiltonian and can be expressed 
in the basis of Slater determinant. It is then convenient to express channel states in (^4)-nucleus 
in the same Slater determinant basis. Consequently, all formal operations involving many-body 
operators and channel/target states become straightforward using the second quantization. 

3.1. GSM-CC equations in the coordinate space 

In order to evaluate Eq. ([2J , we separate the Hamiltonian H = T + U core + V res into basis and 
residual parts: 

H = T + U basis + (V res - U ) , (5) 

where Ub as i s is the optimal potential of ^4-particle system and Uq = Ub as i s — U core . U core is 
the potential generated by the core. The advantage of this decomposition is that V res — Uq is 
finite-range and T + Ub as is is diagonal in the basis of Slater determinants used. 

Let us consider first: |£j ji r») = \£f jf Tf) = \£jr) in Eq. ([2]), as H in this case 
contains infinite-range components leading to Dirac delta's which have to be calculated 
analytically. In order to derive these expressions, we suppose that only a finite number of Slater 
determinants appear in target many-body states. This assumption is always valid because GSM 
eigenvectors expansion coefficients decrease exponentially with the energy of basis scattering 
Slater determinants and, therefore, can be approximated with any arbitrary precision by a finite 
expansion of Slater determinants. As a consequence, the antisymmetry between \u n £jr) and 
|^)^ in Eq. no longer plays any role for n larger than a given n max . This implies that the 
creation operators in Eq. can be replaced by the tensor products. 

It is convenient to have matrix elements {aj3\V res — Uq\j5) vanished when n a > n max (same 
for /3,7 or 8). This is always the case in GSM calculations as one uses a finite model space. It 



is thus convenient to rewrite the Hamiltonian H of Eq. ([5]) introducing an operator acting only 
on the target: 

H = T + U basis + (V res - Uof- 1 + [(V res - Uq) - (V res - Uq)*- 1 } , (6) 

where one defines (V res — Uq)^ 1 by its action on the non-antisymmetrized ^4-body states 
|*)m ® \u n ljmT): 

(V res ~ ^o) A-1 (|*>Af ® K^mr}) = \{V res - U )\^) J M ] ® \u n ljmr) , (7) 

i.e., {Vres— Uq) a ~ 1 is the finite-range part of V re s— Uq acting on (A— l)-body states only. Inserting 
©, ((TJ) in Q, one can show that it is only the sum involving T + C/{, as j s + (V res — Uq) a ~^ and 
n > which generates Dirac delta's in the matrix element (J2J). Consequently, one can rewrite 
the matrix element as a sum of two terms, one which is finite and the other which is infinite: 

(A{{9 f \ J t (rf£jT\}ft A \H\A{m J ' <8> Ir^r)}^) 

Umax 

- £ i^fa) ii„(r/) ({<*/| J ' 5$ (n^jr|}^|T + C/ 6asis + (F res - t7 ) A_1 |{|^) Ji ® Mjr)}j 
+ E u n (ri) ^(r,) <{<*/| J / 8) (u n ^r|}^jT + [/ feasis + (V res - Ub) A_1 |{|*i> Jl ® Mjr)}& A ) 

n>0 



The term in between square bracket is finite and can be calculated using Slater determinant 
expansion of considered many-body states and employing standard SM formulas. As l^j)^. and 

1^ /)m are e ig envec t° rs °f H, the second term of Eq. ([8]) does not vanish if: 

m J A = = m J M ■ 

If \£{ ji Ti) j^z \lj jj Tf), one can show using the orthogonality of \£{ ji Tj) and \£f jf r/) that 
only the sum over rii,nf in Eq. ([S]) is nonzero and can be calculated straightforwardly from the 

Slater determinant expansions of l^i}^. and 
4. Resolution of the CC equations 

Let us consider the scattering ^4-body state decomposed in reaction channels: 



r+oo 
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where c is the reaction channel defined by the (A — l)-body state \^ c )m c an d the one-body 
quantum numbers (£ c ,jc, T c)- u c (r) is the radial amplitude of the c channel to be determined. 
The CC equations in this basis follow from the Schrodinger equation H\$>) = E\$>): 

r+oo 

E / (A{{V e \ Je ®{r ic3cr\} jA \H -E\A{\^ d ) J ^\r'l cl j cl T cl )} jA )u cl {r'){r') 2 dr' = . (10) 



Eq. (|10p is a generalized eigenvalue problem. Indeed, different channels are mutually non- 
orthogonal because target and projectile are antisymmetrized. In order to deal with a standard 
eigenvalue problem, one introduces the channel overlap matrix: 

0(n,c,n',c') = (A{(V c \ J c <g> (n I c j c t\} Ja \A{\* & ) J °' ® \n'l d j d T d )Y A ) . (11) 

Eq. (fit)]) can then be written in a matrix form: HU — EOJJ , where U — {w c (r)} c is 
the vector of considered channels. Introducing: W = O^U, and the modified Hamiltonian: 
H m = 0~zHO~2, one obtains the standard eigenvalue problem: H m W = EW . 

The overlap matrix O is defined in the Berggren basis and can be calculated using the 
Slater determinant expansion of channels. As the antisymmetry acts locally, it is convenient 
to introduce the harmonic oscillator expansion for the finite-range part of 0~2. For this, O is 
expanded in the basis of harmonic oscillator channels to obtain the finite-range part of 0~^: 
A = 0~2 — I d . Then, H m can be separated into long- and short-range parts: 

H m = H + HA + AH + AHA , 

where all terms involving A are expanded in the harmonic oscillator basis. Thus, the added part 
of H in H m can be treated similarly to the short-range residual interaction. 

Using results of Sec. ()3.ip and the transformation described above, one can write Eq. ()10p 
as a system of non-local differential equation with respectively local (V^° c \r)) and non-local 
(vj™ on loc \r,r') Vc) optical potentials: 

g (-<(r) + 4( % + V (r)) + £ {v% c \r)w d {r) + [ (r, r')w d (r') d/j = 

= (E - E Tc ) w c (r) , (12) 

where [i is the reduced mass of the particle and Et c is the energy of \^ c )m c - This set of non- 
local differential equations is then solved numerically using the modified equivalent potential 
method which yields equations local without singularities in the potentials [40J. The initial 
vector U = {u c (r)} c of channel functions is obtained by multiplying solutions of Eq. (I12p : 
W = {w c (r)} c , by O-S: U = W + AW. 

5. Discussion of results for 6 He+p reaction 

Results presented in this section correspond to 6 He target nucleus in 2 states: J w = 0^ and 2^. 
The configuration space for neutrons correspond to 0p 3 / 2 resonance and 17 states of a discretized 
P3/2 contour in the complex k plane. The configuration space for protons includes Opi/2 and 
Ops/2 resonances and 17 states of a discretized contour for each resonance. Moreover, for protons 
we include partial waves s, d, f, g, h which are decomposed using a real-energy contour as in this 
case the resonant poles are high in energy and very broad. 

The s.p. basis in 6 He and 7 Li is generated by the same GHF potential, produced by the WS 
potential of the core and the MSGI two-body interaction between valence nucleons. Parameters 
of the core potential: the radius R = 1.993 fm (1.954 fm), the depth of the central part 
Vq = 50.188 MeV (49.268 MeV), the diffuseness d = 0.497 fm (0.657 fm), and the spin-orbit 
strength V so = 7.488 MeV (8.23 MeV) for protons (neutrons) have been chosen to fit pi/ 2 and 
s i/2 phase shifts in 4 He+p and 4 He+n, and energies/widths of 3/27/ and 1/27/ resonances in 5 He 
and 5 Li. The radius of the Coulomb potential in this calculation is Rc = 1.952 fm. Parameters 
of the MSGI interaction have been chosen to reproduce energies of the states in 6 He (O^ and 
2^) and in 7 Li (3/27/, 1/2]/, 7/2]/ and 5/27/) relative to the 4 He core. Eq. (fT2|) are then solved 
using the multidimensional variant of one-dimensional iterative procedure [40J where the input 
functions for the first iteration come from the diagonalization of H m in the Berggren basis. 
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Figure 2. Relative energy difference between GSM and GSM-CC calculations for He in 
J 71 " = o + (the solid line) and 2 + (the dashed line) states as a function of the depth of the 
WS potential generated by 4 He core. The model space is spanned by 0p 3 / 2 resonance and 30 
P3/2 states of the discretized continuum. As a two-body interaction, the MSGI interaction with 
coupling strengths: V{J = 0,T = 1) = -20MeVxfm 3 and V(J = 2,T = 1) = -15MeVxfm 3 is 
used. 
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Figure 3. The real part of 3/2~ scattering wave function at the CM energy Ecm=5 MeV in 
the entrance channel [ 6 He(0 1 ") (8)^3/2] is plotted for different discretization densities of the p 3 / 2 
contour in the complex fc-plane. 



5.1. Numerical tests 

Before showing results for 6 He+p scattering, it is instructive to compare the spectra of GSM- 
CC calculations with those obtained in GSM by a direct diagonalization in Berggren basis. One 
should remember that the configuration space in both calculations is in general not the same. 
Fig. [2] compares results of the GSM calculation for 6 He with the GSM-CC results for a system 
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Figure 4. Elastic (left panel) and inelastic (right panel) GSM-CC angular cross sections for 
the reaction 6 He+p at the CM energy Eqm = 35.1 MeV. The experimental data are taken from 
Ref. 

5 He(3/2J~)+n in a model space including 0p 3 / 2 s.p. resonance and p 3 / 2 continuum discretized 
with 30 points. One can see that if the continuum coupling is weak, i.e. for deep WS potentials, 
then the relative difference between GSM and GSM-CC calculations tends to zero. On the 
contrary, for weakly bound systems the GSM-CC approach with a limited number of reaction 
channels can be a poor approximation of the GSM calculation in the complete many-body space. 

The discretization density of the Berggren basis is an essential ingredient of the GSM-CC 
calculations. Fig. [3]shows a convergence of the 3/2~ scattering wave function in 7 Li as a function 
of the number of states on the discretized contour. One can see that the scattering wave 
function is fully converged with 17 continuum states. 

5.2. Elastic and inelastic e He+p cross section 

Angular cross-sections calculated in GSM-CC for 6 He+p reaction at Eqm = 35.1 MeV are in 
the qualitative agreement with experimental data[4T] (see Fig. [3]). Quantitative discrepancies 
between theory and experiment may have several origins. Firstly, the excitation energy in Fig. 
[3] is such that the contribution of He core excitations in this reaction cannot be neglected. 
Secondly, the chosen effective interaction is rather schematic. Thirdly and most importantly, 
the target is a weakly bound halo state whereas the reaction product ( 7 Li) in its ground state is 
well bound with strong 3 H cluster correlations due to the proximity of 3 H decay threshold. This 
implies that the GSM-CC configuration space generated by adding one proton to 6 He in 2 discrete 
states (0^~ ground state and 2f resonance) only, may by insufficient to produce configurations 
obtained by a direct diagonalization of 7 Li in the Berggren basis. Indeed, the configuration spaces 
in GSM-CC and GSM become equivalent if all discrete and continuum states of 6 He are used 
to generate 7 Li configurations. Including only few states of 6 He in GSM-CC approach requires 
a strong renormalization of the effective nucleon-nucleon interaction (the optical potentials) 
to compensate for the missing configurations. By construction, such a renormalization will 
generate identical energy spectra for both full space and restricted space GSM-CC calculations. 
However, the reaction cross-sections in these two settings will be different. The comparison 
between energy spectra obtained for the same Hamiltonian either by a direct diagonalization in 
Berggren basis or by solving GSM-CC equations for a limited number of many-body target states 
gives indication of how strong should be the renormalization of effective two-body interaction 
(the optical potentials) due to the neglected channels. In Table [H we compare energy spectra 



Table 1. GSM and GSM-CC spectra for low- lying states of Li are compared with the 
experimental data. For the details of these calculations, see the discussion in the text. 



J 77 


E GSM (MeV) 


Egsm-cc (MeV) 


E exp (MeV) 


3/2~ 


-17.83 


-10.946 


-10.949 


1/2- 


-21.18 


-10.469 


-10.471 


7/2" 


-21.01 


-6.307 


-6.297 


5/2" 


-31.26 


-4.509 


-4.345 



obtained in GSM and GSM-CC approaches using the same MSGI interaction which was fitted 
in the GSM-CC approach to reproduce the low-lying states of 6 He and Li. One may see that 
Li configurations involving 6 He scattering states which are missing in the GSM-CC calculation, 
lead to a dramatic lowering of the GSM eigenenergies and, hence, are essential to understand 
the dynamics in the 6 He+p reaction. 

Results shown in Table [T] demonstrate that none of the CC approaches which neglects 
continuum states of the target can provide a realistic description of the proton scattering on 
weakly-bound neutron-rich nuclei. Similar conclusion about the importance of high-lying states 
in the continuum of a target (A — l)-nucleus on low-lying properties of the (A)-nucleus have 
been made recently in the Shell Model Embedded in the Continuum [42] . 

6. Conclusions 

We have proposed the unified description of nuclear structure and nuclear reactions in the 
framework of GSM which includes continuum couplings in the many-body framework. This 
CC formulation of the GSM can be employed further in the No-Core GSM framework[27j with 
the bare interaction between free nucleons. The great advantage of the GSM-CC formalism is 
that the approximation of neglecting high-lying target states can be checked by comparing the 
calculated GSM-CC energy spectrum with those obtained in (complete) GSM calculation. In this 
way, one can quantify the role of these neglected configurations in the (^4)-particle wave function 
and estimate the scale of renormalization corrections in the microscopic optical potentials. 

In spite of recent developments in ab initio description of nuclear states and progress in open 
quantum system formulation of the nuclear many-body problem, the microscopic description of 
nuclear reactions with weakly bound targets continues to be a formidable challenge as it requires 
including large number of target states to obtain convergent results. In this respect, the task of 
a unified microscopic description of nuclear structure and nuclear reactions with weakly bound 
exotic nuclei remains still a distant perspective. 
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